Spin tunneling in the Kagome antiferromagnet 
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The collective tunneling of a small cluster of spins between two degenerate ground 
state configurations of the Kagome-lattice quantum Heisenberg antiferromagnet is 
studied. The cluster consists of the six spins on a hexagon of the lattice. The 
resulting tunnel splitting energy A is calculated in detail, including the prefactor 
to the exponential exp(— S /fi). This is done by setting up a coherent spin state 
path integral in imaginary time and evaluating it by the method of steepest de- 
scent. The hexagon tunneling problem is mapped onto a much simpler tunneling 
problem, involving only one collective degree of freedom, which can be treated by 
known methods. It is found that for half-odd-integer spins, the tunneling amplitude 
■ and the tunnel splitting energy are exactly zero, because of destructive interference 

ON ■ between symmetry-related (+)-instanton and (— )-instanton tunneling paths. This 

destructive interference is shown to occur also for certain larger loops of spins on the 
Kagome lattice. For small, integer spins, our results suggest that tunneling strongly 

competes with in-plane order-from-disorder selection effects; it constitutes a disor- 
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dering mechanism that might drive the system into a partially disordered ground 
state, related to a spin nematic. 
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I. INTRODUCTION 



Consider a Heisenberg antiferromagnet (AFM) with non-random but competing exchange 
interactions. The classical ground state is often non-trivially non-unique, in having a con- 
tinuous manifold of inequivalent (but degenerate) ground states. However, if one takes 
account of quantum and thermal fluctuations around the classical ground states, the non- 
trivial degeneracy may be broken. The effects of fluctuations can generally be represented 
by an effective "selection" Hamiltonian, which is a function of the classical spin directions 
and "selects" certain ground states (sometimes having long-range order) in favor of others. 
Since long-range order can thus be induced out of an apparently disordered manifold of 
ground states, such selection effects are called "ordering due to disorder" [?]. 

An effect that competes with "order from disorder" selection effects is tunneling between 
different ground state configurations. Tunneling tends to drive the system into a superpo- 
sition of degenerate states, rather than selecting a particular one. Hence, in the regions of 
parameter space in which tunneling events are important, they could suppress "order from 
disorder" selection effects. 

In this paper we study spin tunneling in a 2D quantum Heisenberg antiferromagnet 
on a Kagome lattice [?,?,?,?]. This is a frustrated spin system with a very large ground 
state degeneracy, in which various selection effects have been investigated [?,?,?,?,?]. We 
study tunneling events that involve the rigid simultaneous rotation of small groups of spins 
( "weathervane modes", defined in section ??, in particular a mode involving only the six 
spins on a Kagome hexagon). We find, rather unexpectedly, that the tunneling amplitudes 
are zero when the spin s is half-odd-integer, but non-zero when s is integer. This is due 
to destructive interference between two topologically distinct tunneling paths connecting 
the same initial and final states. The interference occurs when the tunneling amplitudes 
have different topological phase factors [?]. Therefore there should be interesting integer vs. 
half-odd-integer s effects for that range of s-values for which tunneling effects are as strong 
or stronger than selection effects: systems with integer s would have a greater tendency to 
be disordered because tunneling suppresses selection effects, whereas systems with half-odd- 
integer s, where tunneling is absent, would tend to be ordered. 

Apart from studying the role of tunneling in the Kagome lattice, we hope that this paper 
will provide an instructive example of a rather non-trivial spin-tunneling calculation. As is 
customary, the tunneling amplitude of interest is calculated by setting up a coherent-spin- 
state path integral in imaginary time and evaluating it by the method of steepest descent 
[?], which is an expansion in powers of 1/s. Our calculation includes a complete evaluation 

of the prefactor to the exponential e~^°^ h (S Q is the classical action), and a discussion of the 
integration ranges of the spin path integral (these ranges are finite originally, but need to be 
extended to infinity to allow an evaluation of the prefactor). Although the the calculation of 
tunneling rates for spin systems has been of interest in various different contexts [?,?,?,?,?], 
we are aware of only one recent paper where such prefactors are calculated explicitly [?]. 
Moreover, we show explicitly how one may give an exact treatment of the simultaneous 
collective motion of all the relevant spin degrees of freedom, by reducing the problem to one 
involving only a single, collective degree of freedom. The reduced problem can be treated 
by methods well-known from studying a particle in a double-well potential. 

The results of our calculations suggest that for small, integer s, the tunneling amplitude 
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is sufficiently large that tunneling can be regarded as a significant disordering mechanism, 
that tends to drive the system into a partially disordered state, related to a spin-nematic 
state [?,?]. 

This paper is organized as follows: In section ?? we review ground state selection effects 
on the Kagome lattice, and describe the hexagon tunneling event that is to be studied in later 
sections. In sections ?? and ?? we study a simple model Lagrangian, chosen such that it will 
be of use for the subsequent study of the Kagome system. Specifically, the calculations of the 
classical action S Q (section ??) and the tunnel splitting energy A (section ??) are presented 
in detail. Section ?? shows how the full hexagon tunneling problem can be mapped onto 
the simple model problem studied in earlier sections, and contains the main result of this 
paper, eq. (??). In section ?? the occurrence of destructive interference during the tunneling 
of larger sets of spins on the Kagome lattice is discussed. There are four appendices. In 
appendix ?? an estimate is made of the size and shape of the coplanarity potential that 
we employ in section ?? to study the hexagon tunneling problem. Appendix ?? contains 
a summary of results that are well-known for a particle tunneling in a double well and are 
needed in section ??. 



II. THE KAGOME ANTIFERROMAGNET 

The Kagome lattice is a 2D triangular lattice with lattice constant a Q , where sites have 
been removed at all sites of a triangular superlattice, with lattice constant 2a Q (see fig. ??). 
The quantum Heisenberg AFM on this lattice has the Hamilton operator 

H = J'£si-s j , (J>0) (2.1) 

where 3j is the spin vector operator for the i-th spin and the sum runs over all near- 
est neighbors. Coherent spin states (see, e.g. [?]) may be used to discuss this system in 
classical terms. Associate the coherent spin state \(pi,0i) with the i-th spin, where the 
spherical coordinates ((fii,9i) = fi« define a unit vector rii (fig. ??a). The well-known prop- 
erty {<t>i,9i\si\<t>i,9i) = srii, allows one to interpret rij as a "classical spin vector". Also, a 
"classical" Hamiltonian H may be introduced as the coherent-spin-state expectation value 
(H): 

H= (n^lj^ini^)) = s 2 J'£n i -n j , (J>0) (2.2) 

When referring to the "energy" of a state, we shall always mean this expectation value TC. 
Likewise, the term "ground state" will not be used to refer to an actual eigenstate of the 
operator H, but to a state that minimizes H. It is this Ti that is used in discussions of the 
classical Kagome antiferromagnet. 

The energy is minimized by any configuration in which the total spin on each elementary 
triangle of the lattice is zero [?,?]. In such states, the spins of any given triangle of the 
Kagome lattice lie in one plane, with relative angles of 120° (see fig. ??b), forming a rigid 
unit in spin space. Because of the many ways of fitting such triangles together, the distinct 
classical ground states form a manifold with a dimensionality proportional to the system 
size. 



3 



A. Selection effects 



It is convenient to describe the classical ground states with reference to a coplanar ground 
state, in which all spins lie in the same "reference plane". All the coplanar ground states 
are constructed by labelling the sites by letters A,B and C, such that each triangle has one 
of each, and then replacing the three kinds of letters by spins in three directions differing 
by 120° angles [?] (fig. ??b); these states correspond one-to-one to the ground states of a 
3-state Potts AFM on a Kagome lattice. 

All non-coplanar classical states can be generated by continuous distortions of coplanar 
states, henceforth to be called "Potts states", without crossing energy barriers [?]. For 
example, any hexagon of six spins that is labelled only by two letters (e.g. ABABAB in 
fig. ??) can be rotated as a rigid unit in spin space by an angle tp around the C-direction 
without any cost in classical energy; this has been called the "weathervane" mode by Ritchey 
et. al. [?]. However, by expanding in spin waves about any given ground state, it has been 
shown that all non-coplanar states have a larger zero-point energy than the coplanar Potts 
states. This selection effect can be characterized by a parameter ~ 0(sJ), and ensures 
that the true ground state will be coplanar [?,?,?,?,?]. 

Further study, using a self-consistent approach, yields a smaller "in-plane" selection 
energy, of order J c ~ 0(ss J), which favors among all possible Potts states a particular 
state with long-range order, the so-called a/3 x a/3 ground state [?,?,?], depicted in fig. ??. 
In this state, every hexagon is labelled by just two letters, e.g. BCBCBC. To encode the 
most important effects of this "in-plane" selection, it is convenient to define a "chirality" 
variable centered on each triangle, equal to —1 or +1 depending on whether the sites take 
values ABC in the clockwise or counterclockwise sense. In the a/3 x a/3 Potts state, the 
chiralities on neighboring elementary triangles are antiferromagnetically ordered. Hence the 
selection energy can be expressed as an effective antiferromagnetic coupling, of strength J c , 
say, between the chiralities. 

B. Competition between selection effects and tunneling 

The above considerations neglected the possibilities of large-scale fluctuations and tun- 
neling between classical ground states. Clearly the smallest object which can tunnel is the 
"weathervane mode": Consider as initial state, the a/3 x a/3 Potts state. If the six spins 
on an ABABAB-hexagon rotate by 180° around the C-direction (take the latter to define 
the z-axis), another Potts state, |/), with a BABABA-hexagon is reached (fig. ??). In the 
absence of in-plane selection (J c = 0), \i) and |/) would be degenerate. Tunneling between 
|«)and \f) would tend to drive the system into a superposition ± |/)), with energy 

E ± A (where A is proportional to the tunneling amplitude). Now, in the presence of 
in-plane selection effects, the energies of \i) and |/) differ by 12 J c , since the chiralities on 
all six triangles bordering the hexagon of \i) are opposite to those of |/). If 12 J c > A, the 
tunneling amplitude is very small, and the system selects \i), the a/3 x a/3 state. If, on the 
other hand, 12J C < A and the tunneling amplitude is appreciable, the system can lower its 
energy by adopting the superposition — |/))- 

Such a hexagon tunneling event can clearly take place starting from any Potts state \i) 
that contains an ABABAB-type hexagon, not just from the a/3 x a/3 state. Thus, tunneling 
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constitutes a disordering factor that competes with in-plane selection. If such hexagon 
tunneling events occur with large probability throughout the Kagome lattice, a ground 
state will result that is at most partially ordered (related to a spin nematic) (see section 
??). 

C. Effective hexagon Hamiltonian 

We shall study the hexagon tunneling event described in the previous section in the 
absence of in-plane selection (i.e. taking J c = 0), but in the presence of a coplanarity barrier 
of size Jb, that tends to keep the spins aligned with the reference plane. Our aim is to 
compute the tunnel splitting energy A. Since Jb ~ sJ, we have Jb <C s 2 J, when s >> 1, 
and consequently we assume that only the six hexagon spins will move significantly during 
the hexagon tunneling event. Hence we take all other spins in the lattice to stay fixed, and 
adopt the following Hamiltonian for the six-spin hexagon system: 

T^hex = 7~^AFM + H cop , (2.3) 

5 

where Hafm = s 2 J^][n r z + n i+1 -z + n r n /+ i + §] , (2.4) 

1=0 

5 

H cop = J b f((Pay) , where y? av = 1/6 ^ ((pi - (initial)) . (2.5) 

1=0 

The index / labels the six spins around the hexagon and is defined modulo 6 (see fig. ??). 
All spins are written in terms of the same coordinate system, in which the z-axis points in 
the C-direction, the xz-plane coincides with the reference plane, and tp is measured from the 
positive x-direction. Hafm contains simply those terms from eq. (??) that involve the six 
hexagon spins and their six C-type nearest neighbors. A constant has been added to ensure 
that Hafm = in the initial and final configurations \i) and |/). 

H cop represents the coplanarity barrier that opposes hexagon tunneling. The function 
f(ip) describes the shape of the barrier; it is of order unity and is sketched schematically 
in fig. ??. Its argument in eq. (??), y? a v, is a measure of the deviation of the plane of the 
near-rigid hexagon from the reference plane (in which (p B (—T/2) = 0, v?a(— T/2) = ir). 

Symmetry about the reference plane ensures that /(</?) has the properties 

f(mir) = , f(ip) = f(—(p) = f((p + rnn) for any integer m . (2.6) 

At this stage it is not necessary to specify f(tp) in more detail. An estimate of the actual 
form of 7i cop is made in appendix ??, and summarized in eq. (??). The suitability of using 
such a barrier is further discussed in section ??. 

Intuitively, the six hexagon spins are expected to rotate collectively, almost as a rigid 
unit, maintaining the mutual coplanarity and relative angles of 120° required to minimize 
Hafm- In particular, for each of them cos#; ~ —\ throughout the tunneling event, which 
is why we have taken Ti, cop to be independent of Q\. Also, the expected mutual coplanarity 
of the tunneling spins is the reason why we have written TC cop as a function only of ip av 
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(and not some more general function of the six ip^s [?]). These assumptions are found to be 
justified in section ??. 

Evidently, due to reflection symmetry in the reference plane, two kinds of tunneling events 
between \i) and |/) are possible. They differ from each other only in the sense of rotation 
about the z-axis (<p <-> —<p), and we shall call them (+)-instantons and (— )-instantons. 

We shall show in section ?? how the hexagon tunneling problem defined above can 
be mapped onto a much simpler model problem. This involves only a single, collective 
spin degree of freedom with effective spin 6s. Its coordinates, to be denoted by ($ ,6 ), 
are formally defined in eqs. (??) and (??) and are suitable averages of the six individual 
(<£>i,0j)'s. The effective Hamiltonian turns out to be (see eq. (??)): 

H eS = 12s 2 J(cos O + \f + J b /($ ) . (2.7) 

Rather than proceeding with the mapping of eq. (??) onto eq. (??) right away, in the 
next two sections we first discuss the simple model problem (based on eq. (??)) in detail, 
to establish notation and introduce the tools needed to calculate the tunnel splitting energy 
A. 

III. A SIMPLE MODEL SYSTEM 

In this section we introduce a simple model problem, which will be used to illustrate how 
a calculation of the tunnel splitting energy can be carried through. It is also the system 
onto which the Kagome tunneling problem that is studied in section ?? can be mapped. We 
set up the relevant path integral, calculate the classical action and discuss the equations of 
motion and the typical form of the tunneling path. 

A. Model Hamiltonian 

The model system is defined by the following Euclidean single-spin Lagrangian: 

C = -ihns(p(z -1) + H (3.1) 
where H = a(z - z g ) 2 + bf(<p) ; (with b < a) , (3.2) 
and z = cos# . (3.3) 

This Lagrangian has been written in terms of the imaginary time parameter t — it (hence 
"Euclidean"), since this is convenient for the calculation of tunneling amplitudes. The 
spherical coordinates (ip, 9) = Q define a unit vector n and label a coherent spin state \ip, 9) 
for a particle with spin ns. The dot on <p means d T (see [?] for a discussion of the origin 
of this term). The integer n is introduced in order to accommodate the possibility of a 
collective degree of freedom with effective spin ns. For the purpose of describing a single 
spin degree of freedom, take n — 1. 

The "classical" Hamiltonian is the expectation value H = ((p,9\H\(p,9) of the quantum 
operator H. Clearly H has been chosen to have the same form as H c g of eq. (??); the 
constants a and b are of order s 2 J and J b respectively (see eq.(??), with b a. The dominant 
term in Ti is an easy-plane anisotropy which would make every angle on the "latitude line" 



6 



z = Zg on the unit sphere be degenerate; it mimics Hafm in eq. (??) (which forced all six 
hexagon spins to have cos#; = — |). In the other term, f((p) is the same function as that 
in eqs. (??) and (??); it introduces a small anisotropy within the degenerate subspace and 
mimics Ti, cop . There are two degenerate ground state configurations, \i) = \2miir,z g ) and 
\f) = |(2m 2 + 1)tt, z g ), with mi and m 2 arbitrary integers (the different values of mi and m 2 
describe the same physical state, of course). Contours of constant H are shown in fig. ??. 

B. Tunneling amplitude 

The tunnel splitting energy A that arises due to tunneling of the spin direction between 
\i) and |/) is proportional [?] to the tunneling amplitude 

U fi = (f\e-^ T / n \i) = J IXle 3 '* , (3.4) 

rp Icy 

where S = J_^ 2 drC is the Euclidean action, VQ is the path integral measure (discussed in 
section ??), and T is a large time. Such an amplitude can be approximately evaluated by 
the method of steepest descent: 

Ufi =Y.M {l) e- S ° l)/h = uf . (3.5) 
i 

Here Sj^ is the action evaluated along the Z-th "tunneling path", which is a solution to the 
Lagrangian equations of motion and for tunneling problems is in general complex. It will 
always be denoted by overlined variables, e.g. (lp( l \z^). The index / allows for the possibility 
of different tunneling paths satisfying boundary conditions that differ in the indices m\ and 
m 2 (but that all describe physically the same initial and final states) [?]. The prefactors 

= JvQe-^-^^ . (3.6) 

are usually evaluated only (if at all) to lowest order in the steepest descent method, by 
transforming to the fluctuation variables (5ip^ l \ 5z^) = (</? — Tp( l \ z — z®) and keeping only 
the lowest term in the expansion S - S® = 5 2 S (l) + 5 3 S (l) .... 

In the present case, all tunneling paths connecting \i) at r = —T/2 to |/) at r = T/2, 
can be constructed from two very simple paths, to be denoted by (Tp ± ,z ± ). The first, for 
which y5 + (r) e [0, n], we call a (+)-instanton, the second, for which Tp~ = —Tp + (r) G [0, — n], 
a (— )-instanton. They differ solely in the sense of (^-rotation and are sketched in fig. ??. All 
other tunneling paths that approximately satisfy the equations of motion and contribute to 
eq. (??) are multi-instanton paths. They consist of n\ (+)-instanton and n 2 (— )-instanton 
events, with n\ + n 2 = odd, all assumed widely separated (relative to their characteristic 
width) in time, and following each other in arbitrary order (this is the so-called dilute-gas 
approximation, see Coleman [?]). In the following, we consider only single-instanton events 
(i.e. I — > ± in eq. (??)); the effects multi-instantons are taken into account in appendix ??. 

The symmetry H((fi, z) = 7i(—(p, z) of the Hamiltonian allows one to conclude immedi- 
ately that = \Ujj\. This is intuitively obvious, and can be proven to hold exactly to 
all orders of the steepest descent approximation (i.e. to all orders in 1/s) (see section ??). 
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Intuitively speaking, the symmetry of 7i ensures that the local neighborhoods of the two 
tunneling paths (lp + ,z + ) and (Tp~,~z~) are identical for the two paths, so that the local 
shapes and sizes of the barriers (which determine Re[«S®]) and the local fluctuations around 
the tunneling paths (which determine \Af^\), are identical. However, the amplitudes 
can differ by a phase, which may lead to destructive interference between them. 



C. Classical action, continuing to complex coordinates 

Let us find , the classical action for single-instanton events. Since energy is conserved 
along any path that extremizes the action, TC(Tp ± ,z ± ) = 0. Solving for z* 1 as a function of 
Tp, one obtains 



iZ ± = z±- Zg = ±ijbfa y/ffi) , (3.7) 

with Z ± real. The Euclidean action is easily evaluated along these paths by changing 
variables from r to Tp ± and using f(^p + ) = f(W~) : 

rT/2 , 

S± = -ihns / rfr(9 T ^ ± ) (z q + iZ - 1) + 

J-T/2 

= ^fihnsn(z g — 1) + hnsyfb/a f dip^J f(lp + ) . (3.8) 

J 

Note that the quantity \Jb/a plays the role of an effective barrier size for the tunneling 
process. This is intuitively plausible. The barrier height in the direction is measured by 
b. When (z — z g ) is real, a measures the steepness of the valley in the z-direction. When 
(z — z g ) becomes purely imaginary ((z — z g ) 2 — > —{z — z g ) 2 ), the ^-valley turns into a ridge 
(see fig. ??). The motion occurs along a constant energy contour around this ridge (see 
fig. ??b), and 1/a measures the width of the ridge. Hence, b/a measures the effective barrier 
size. 



D. Possible cancellation due to phase factors 

We show in section ?? that J\f + = N~ (under some technical assumptions, explained 
there). The sum of the two amplitudes is therefore simply 

U% + U£ = U%{1 + e S °- S °» n ) 

= t^(l + e- ftri ), (3.9) 

where the constant / = 2ns(z g — 1) is not necessarily an integer. The relative phase of e~ ml 
between the amplitudes of two tunneling paths that connect the same initial and final states 
has a well-known geometrical interpretation. It is related to the area enclosed on the unit 
sphere between the two paths ((^ + (r), Re[^ + (r)]) and ((/2~(r), Ke[8~(r)]). (See [?] for more 
comments on this aspect). Now, if I is an odd integer, the (+)-instanton and (— )-instanton 
paths interfere destructively and their amplitudes add to zero: + Uj { = [?]. Note 
that as long as the barrier does not violate the symmetry between (+)-instantons and (— )- 
instantons, this cancellation does not depend on the particular shape of the barrier, since 
the shape function f(<p) only affects the real part of <S G , which cancels out in — S~. 



8 



E. Equations of motion and tunneling path 

It was not necessary to solve the equations of motion to obtain the explicit expression 
eq. (??) for the classical action. For future reference, the equations of motion are: 

-ihnsit = bf'(Jp) , (3.10) 

— ihnsTp + 2aiZ = 0. (3-11) 

The prime on / means 8/ dip. Note how the % of %Z and that originating from t — it combine 
to result in consistent equations of motion (this is a direct illustration of why it is useful to 
employ imaginary time in tunneling calculations). Eliminating Z, we get 

if = ^f'm, where t? = ^ . (3.12) 

The constant d has the dimensions of inverse time, and 1/d characterizes the width of an 
instanton. To illustrate the nature of the solutions of eq.(??), consider the simple case where 
f(ip) = sin 2 if. Then eq.(??) is just the sine-Gordon equation, whose solutions, in the limit 
T — > oo, are: 

JP + {t) = 2 arctan(e d *) , TpT( T ) = -p + (r) = p + (-r) - tt . (3.13) 

These represent (+)-instanton and (— )-instanton events and are sketched in fig. ??. The 
shape of an instanton changes quantitatively, but not qualitatively, if a different barrier 
shape function is used. 



IV. THE TUNNEL SPLITTING ENERGY A 

For integer s, for which the (+)-instanton and (— )-instanton amplitudes have the same 
sign and interfere constructively, the resulting tunnel splitting energy A can be calculated 
from the prefactor of eq. (??): 

M ± = Jvtle-( S - S± ° )lh . (4.1) 
This can be accomplished via three steps: 

(i) Extending the integration ranges for ip and z in the path integral to infinity; 

(ii) Integrating out the z-degree of freedom; 

(ii) Using standard methods to evaluate the resulting 93-path integral. 



A. Extending the integration ranges 

The path integral expression eq. (??) for the tunneling amplitude (in real time) can be 
arrived at by the usual procedure of discretizing time and inserting completeness relations 
in spin space at each time slice (see, e.g. [?,?])• This procedure leads to a formal expression 
for the measure (appropriate for any spin problem): 
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r /2 s 4- 1 \ N N r 2v r 1 

j Vn(t) = lim^ [-^-) II I M^) J dz(Tj) , (4.2) 

where e(N + 1) = T. What distinguishes this measure from the ones usually encountered 
in particle tunneling problems is the fact that <p and z have finite integration ranges. In a 
somewhat cavalier fashion, we extend these to [— oo, oo] and [— oo, oo] and absorb the change 
in normalization by multiplying by an extra overall normalization factor C. The extension 
of the range for tp from [0, 27r] to [— oo, oo] is very natural - it allows for motion in which 
the spin direction rotates around in the same direction many times. The justification for 
extending the integration range [— r, r] for z from r = 1 to r = oo is more tenuous. The 
core of the argument is the following assertion, proven in appendix ??: In the presence of 
a 2 2 -term in the Hamiltonian, the value of r determines the degree of non-differentiability 
( "wildness" ) of the </?-paths that result after z has been integrated out: if r = oo, the <£>-paths 
are Brownian motion paths; if r = 1, they are much more ill-behaved than Brownian motion 
paths. Therefore, changing the integration range for z from r = 1 to r = oo is equivalent to 
restricting attention to Brownian motion y?-paths instead of a larger class of paths that are 
much more ill-behaved. It is argued that this should not have a noticeable effect, for the 
following two reasons: if one is interested mainly in the effect of small fluctuations around 
the classical path (as, for example, when calculating the tunnel splitting energy), physically, 
one only expects some smooth paths in the immediate neighborhood of the classical path 
to be important. Also, the additional paths that are formally included in the path integral 
when the integration ranges are extended from r = 1 to r = oo are far from the tunneling 
path and therefore only make an exponentially small contribution to the path integral. 

In the course of making a change of variables (ip, z) — > (Sip, Sz) in the path integral in 
order to evaluate J\f ± , the overall constant C will be multiplied by a Jacobian factor J N , 
which may be infinite in the limit N — > oo. To obtain finite answers, we stipulate, as is usual, 
that C be chosen such that the final path integral for 5<p (after the ^-dependence has been 
integrated out), should have the same normalization as that employed by Coleman in his 
discussion of instantons [?], the results of which we intend to use. This often-used procedure 
may seem somewhat arbitrary, but the fact that the path integral is defined as an infinite 
product of integrals, each of which may produce finite prefactors when being manipulated, 
leaves one no choice but to absorb all infinities in a single appropriately chosen constant C . 
This much having been said, we henceforth pay no attention to normalization constants or 
Jacobians. 



B. The relative phase of the prefactors N + and N 

Before integrating out z, let us investigate the relative phase of M + and Af~, the (+)- 
instanton and (— )-instanton prefactors. The quantity S — S Q that appears in the integrands 
of the prefactors in eq. (??) can be written in terms of the fluctuations around the tunneling 
path, (Sip, 5z) = (<p — Jp, z — z), as 



S-S 



dr 



-ihnsSipSz + a5z 2 + b ^ 



n=2 



gn 

dip n 



5<p r 



(4.3) 
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Evidently, because f(tp) = /(—</?), the following symmetry relations hold (the square brack- 
ets denote a functional dependence, "*" denotes complex conjugation): 



(S - S )[<F, S<p, Sz] = (S- S )*[ip + , Sp, Sz] , (4.4) 
(S - S )[Tp~, -6<p, -Sz] = (S- S )[Jp + , 6<p, Sz] . (4.5) 

Now, after the change of variables (p, z) — > (Sip, 8z) in the path integral (??), Sp and Sz are 
dummy variables that are integrated over. Since we extended the integration ranges for <p 
and z to [—00, 00], the integration ranges for Sip and Sz are symmetric around Sp = and 
Sz = 0. Hence the following conclusions follow immediately: 

(i) Relation (??) implies that M~ = M + * . 

(ii) Relation (??) implies that Af~ = Af + . 

As discussed in section ??, the extension of the integration range to [—00, 00] is on 
somewhat less firm ground for z than for ip. Hence conclusion (ii), which holds only if the 
^^-integration range is symmetric about Sz = 0, is in a sense "weaker" than conclusion (i). 
However, even if the original ^-integration range of [—1, 1] is retained, the error in the relation 
A/" + = M~ is expected to be exponentially small. For example, for the which 
we shall use in section ??, only paths for which Sz G [\, |] break the Sz <-> —Sz symmetry, 
and these paths deviate so strongly from the tunneling path (for which Re[z] = z g ) that 
their contribution to the path integral is exponentially small. (Of course, this argument 
breaks down when z g = ±1.) 

Finally, note that within the lowest order of the steepest descent approximation, in 
which one keeps only the n = 2 term of the infinite sum in eq. (??) (to be indicated by a 
superscript (2) in eq. (??) below), the relation Af + = Af~ holds, independent of the range 
of ^-integration. The reason is simply that since 

(S-S ) {2) = Jdr[-ihsSpSz + aSz 2 + \bf"{p)Sp 2 ] , (4.6) 

where /" = d 2 f /dp 2 , we have 

(S-S )^\p + ,Sp,Sz} = (S-S f\p-,Sp,Sz] . (4.7) 



C. Calculating A 

We henceforth restrict attention to the lowest order of the steepest descent approxima- 
tion, in which A/" + = A/"~ = M . It is straightforward, starting from eq. (??), to perform the 
Gaussian integral over Sz to arrive at 

AToc Jv5pe-J dT52C * /h , (4.8) 
where S 2 C V = \bf\Tp){Spf + ^^(^) 2 • (4.9) 
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Now note that the 8 2 C lf of eq. (??) is also the second variation of the following effective 
Euclidean Lagrangian for ip: 

C, = ±^<p 2 + bf(<p). (4.10) 

The Lagrangian equation of motion resulting from £ v is just eq. (??), and hence the S Q 
corresponding to C v will be equal to the Re[«S ] found earlier. Consequently, both Af and 
Re[«S ] for the original system are equal to those arising from C v . It follows that A, too, 
can be calculated directly using C v . Furthermore, C v is quadratic in Sip and therefore the 
methods discussed by Coleman [?] can be used to calculate the tunnel splitting energy. 
Coleman's methods, which are summarized in appendix ??, readily lead to the following 
expression for A (compare eqs. (??) and (??)): 

A = 2v?as(7ms)-5(6/"(0))i(2a)3 (e~ S ° /h + e~ S ~° ,h ^j . (4.11) 

Here the constant y2 as is to be read off from the asymptotic behavior of the tunneling paths, 
which, as Coleman shows, can always be written in the form (compare eq. (??): 



J2abf"(0) 

lp ± (r)^±TiT^e- u>T } withcc^L— , (4.12) 

Eq. (??) is the main result of this section. We emphasize once again that A is strictly 
zero if 5+ — S~ = ihlir, with I an odd integer, since then (+)-instanton and (— )-instanton 
events interfere destructively. 

To find (pus explicitly, one needs some knowledge of the asymptotic behavior of the shape 
function f(ip) of the potential as (p — > 0. According to eqs. (??), (??) and (??) in appendix 
D, the form of f(ip) that is applicable to the hexagon tunneling problem to be studied in 
the latter parts of this paper is: 



/(¥>)=/ V^ff + V* 2 " /(</?*) , = V mod tt G (-vr/2, tt/2] , (4.13) 



where /(</?) = | sin<^|(l + |sin 2 (/9)2 — J sin 2 (p . (4-14) 

Here J ~ 0.42 and <p* = 0.14 s -1 / 3 . In the limit s ^> 1, <p* can be treated as a small parameter 
which characterizes the curvature of the potential at <p — 0, since /"(0) = (l/<p* — 2J). 
"Energy conservation" along the tunneling path implies for the C v of eq. (??) that 

|^V-^)=0. (4.15) 
Integrating this equation gives, for an (— )-instanton, 

TO(m -i = -/7**, (4.16) 

J-n/2 [2f(<f)}2 
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where — it + Tp T is the angle reached at time r if Tp = at r = 0. This equation is to be 
solved for ~tp T as a function of r, in the limit Tp T — > 0. In particular, we take ^ r <C r^y?*, where 
77 <C 1 is an arbitrary small parameter (e.g. r\ = 0.01). 

It is convenient to split the integral into two parts by writing 

M/ ff (0))-* = *W2) - , (4-17) 

where F(y?) = T — , (4.18) 

and the property f(ip) = f{tp + n) has been exploited. Since Tp T <C rjtp* <C 1, one may 
evaluate F(tp T ) analytically by using the asymptotic form of the integrand, namely f{tp — > 
0) ~ (p 2 (l/2(p* - J), with the result 

F&r) = ( iXjtp* ) 2 (lllV? " ■ (419) 

Using this result in eq. (??) and solving for Tp T gives 

tp T = tp^e-™ , (4.20) 
1 -2J<p*\ 



where tp^ = exp 



C 



tp* 



(4.21) 



and C = F(tt/2) + ( J * ln^y?* (4.22) 

\1 — 2Jtp* J 

The constant C is independent of r], since the ^-dependence of F(tt/2) cancels that of 
the second term. Moreover, C is only very weakly dependent on tp* and hence on s, and 
approaches a constant value as s — > oo. (Numerically it is found that C changes from 1.1 
to 2.1 as s changes from 1 to oo.) Hence, recalling that tp* = 0.14s~3, we conclude that in 
the limit s — > oo, 

V?as^e 5 - 6 ^ . (4.23) 

This expression diverges as s — > oo, but in eq. (??) it premultiplies an expression that tends 
to zero even faster as s — > oo (in section ??, we find S Q ~ s 1//2 , see eq. (??) and eq. (??)). 

V. SPIN TUNNELING IN THE KAGOME LATTICE 

The stage has now been set for the study of the hexagon tunneling event on the Kagome- 
lattice, for which we adopted the Hamiltonian Hhex defined in eqs. (??) to (??). Rewriting 
TCafm in terms of tpi and z\ = cos#;, the Euclidean Lagrangian to be studied is 

5 

£ = -lhstpi(zi - 1) + Jbf(<fa.v) (5.1) 
1=0 

5 

+ s2 J Yl\. z i + + z i z i+^ + C i ^(fi - Vi+i) + |] > 

Z=0 
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where C, = sin 0, sin 0,+i = {(1 - zf)(l - zf +l )}^ . (5.2) 

The tunneling event was described in intuitive terms in section ??. The initial and final 
states are: 

(pi(-T/2) = (0,7r,0,7r,0,7r) ; zi(-T/2) = z g . (5.3) 
|/>: ^(T/2) = (7r,0,7r,0,7r,0); ^(T/2) = z g . (5.4) 

Here z 9 = cos27r/3 = — \. The motion of each of the six spins will be roughly analogous to 
the single-spin motion described in section ??. 

A. Collective coordinates 

We now introduce a coordinate transformation to a set of collective coordinates, (<$>i, Zi), 
which ultimately allows us to map the hexagon tunneling problem exactly onto the simple 
model problem of section ??, eqs. (??) and (??): 

z l -z l {-T/2) = j^e Mk "iZ k , (5.5) 

fc=0 

<Pi-<Pi(-T/2) = j2e~ Mk/3 *k, (5-6) 

fc=0 

where = , Z\ = Z_ x . (5.7) 

The conditions (??) ensure that (fi — (fi(— T/2) is purely real and zi — zi(—T/2) purely 
imaginary. The reason for the latter requirement is the same as that encountered in the 
discussion of the model tunneling problem (section ??): it would otherwise be impossible to 
satisfy 7i — 0. Making zi — Zi(—T/2) imaginary flips the sign of the z 2 -terms, thus turning 
the z-valley into a ridge around which the tunneling path can proceed along a constant 
energy contour. 

The new coordinates are essentially discrete Fourier transforms of the original ones with 
respect to I. (Indices such as I and k will be used interchangeably both as "position" labels 
(for if and z) and as Fourier labels (for $ and Z).) The zeroth ones, ($ D , Z a ), are just 
the average values of the old coordinates and can be thought of as the coordinates of a 
collective degree of freedom (of spin 6s, as it turns out). They will be referred to as the 
tunneling coordinates, since $ G tunnels from at r = —T/2 to ±7r at r = T/2, and iZ Q is its 
conjugate momentum. The other coordinates ($/, Zi), I — 1, ... ,5, will be called transverse 
coordinates, since they will be shown to be strictly zero along the tunneling path, and hence 
describe fluctuations that are orthogonal to the tunneling path. 

In order to calculate the tunnel splitting energy A, one needs the classical action S ± 
and the prefactor J\f ± . In the following two sections, we show that both these quantities can 
be obtained from a rather simple effective Lagrangian, C Q , that depends only on the tunnel- 
ing coordinates (Q ,Z ), so that the problem is substantially simplified. This Lagrangian, 
defined in eq.(??) and given explicitly in eq.(??), turns out to have just the form of the 
Lagrangian introduced in section ??, with (<3> , Z Q ) corresponding to if, z) there. 
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B. Classical action (Kagome lattice) 



The classical action is completely determined by the tunneling path, for which, as in 
section ??, we again use overlined variables, (ifi,zi) or (<I>j,Zj). This path is determined 
by the equations of motion, which, by the chain rule, can be written in terms of the new 
coordinates as: 



dC\ ( d(fi 



dipi) t \d$ k J t 
K dz l )\dZ k ) t 



'dC_ 

d<pi 



' dipt 



0, 



(5.8) 
(5.9) 



with boundary conditions $(— T/2) = and Z(—T/2) = 0. The bracket ( ) t indicates that 
the derivatives are to be evaluated along the tunneling path. These equations can be solved 
by the following "tunneling path Ansatz" : 



($,,Z,) = (0,0), for/ = l,...,5. 
To see that this Ansatz works, consider eqs. (??) and (??) term by term: 



(5.10) 



(5.11) 



'd£\ ((hpA 



[s 2 J [-C, sin(y>, - ip l+1 ) + Ci-i sin(^_ 1 - ip t ) 
+ |^/ , (^av)}e-^ fc / 3 ; 



(5.12) 



dzi) t \dZ k ) t 



s 2 J 



— ihs(p l > ie 1 



T/fc/3 



(5.13) 



with 



' dC m \ 
> dz i It 



Zm{l — Z m+1 )8 m i + Z m+ i(l — Z m )5 m+ n 



(5.14) 



Now express (ipi,zi), wherever they occur, in terms of (<&i,Zi). Using the tunneling path 
Ansatz (??) repeatedly, it readily follows that 



zi = iZ Q + z g = Z c ; ip t = $ c + (p (-T/2) ; 



(5.15) 



whence 



x""-^ / dC m \ 

£^ V dz i it 



<Pl - <p l+1 = 7T 



= -2Z c . 



(5.16) 
(5.17) 
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This means that eqs. (??) to (??) depend on the index I only through the factors of e ±lwlk ^ 3 . 
Since Yh=o e ±l7vlk ^ = 6S k o, these equations simplify considerably to become 

- ihsiZ~ 5 Qk = \ J b f'(§o) S ok , (5.18) 

ihst o 5 0k = As 2 JiZ 5 Qk . (5.19) 

where z g = — \ has been used. Evidently, the equations of motion for the transverse coor- 
dinates (k 7^ 0) trivially reduce to zero, and only the tunneling coordinates $ and Z Q have 
non-trivial behavior. 

This allows us to introduce a considerable simplification: Write £ for the function that 
results when C is expressed in terms of the new coordinates: 

Z ; $ 5 , Z 5 ) = C ({¥>*({$,}), ZkdZt})) . (5.20) 

Define an effective Lagrangian C , that depends only on the tunneling coordinates, by the 
relation 

£ ($ , Z ) ee Z o] 0, 0; . . . , 0; 0) . (5.21) 

Now, the behavior of the tunneling coordinates $ and Z Q can be found directly via the 
effective Lagrangian C Q) instead of the full C. This follows because the equations of motion 

for C are exactly eqs. (??) and (??), essentially by definition, since , evaluated 

mdmg t0 *. — g Path Ansatz e, (TT). - e qual t „ (^)^ ^ Similar 

comments holding for (^ d ^ d £ and (^j^J ■ Furthermore, since S Q depends solely on the 

tunneling path, it too can be calculated directly from C Q . The effective Lagrangian C Q can 
be found from eq. (??) to be: 

£ ($ , Z ) = -ih6s(6 )(Z c - 1) + 12s 2 J{Z C - z g ) 2 + J b f{^ ) . (5.22) 

It is from this expression that H c s of eq. (??) was obtained. Evidently C has just the form 
of the Lagrangian of eq. (??), with 

n = 6; a = 12s 2 J ; b = J b ; z g = -\. (5.23) 

Consequently, all the results from section ?? are applicable. The classical action, from 
eq. (??), is 



St = ±9ihs7r + h S] j j± J o n d^f(ip+) . (5.24) 

From the above expression for the classical action one can immediately read off the most 
striking result of this paper: 



e -S + jh + e S~/n = J -2exp [-^/3J b /J j^d^ f(^ + )\ if s = integer , ^ ^ 

) if s = half-odd-integer . 

Thus, if s is half-odd-integer, the (+)-instanton and (— )-instanton amplitudes interfere 
destructively, and the total tunneling amplitude is zero (for a discussion and extension of 
this result, see section ??). 
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C. The prefactor (Kagome lattice) 



The calculation of the prefactor involves the evaluation of the path integral 

5 

JflV^e-f^ 52 ^ . (5.26) 



=o 

where the second variation of the action is 



(5.27) 



The deviations around the tunneling path, defined by 



b Zl = Zl -z x = e Mk / 3 5Z k , (5.28) 

Sw^w-n ee e—' fc / 3 5$ fc , (5.29) 

are taken to be purely real (it is only the tunneling path itself that has to become complex 
to minimize the action). Thus, their Fourier transforms obey the conditions 5^ = <5$_z 
and 5Z* = 5Z_i. 

Now, it can be verified that 

( d 2 C \ ( d 2 C \ 

-ihsd mn ; — = , (5.30) 



dip m dz n ) t \dip m dz n/ 
"~ L " ' -\Js 2 {25 mn ^^ + {5 mn+l + 5 m+ln ){l-^^)\ (5.31) 



d 2 c 



d(p m dip. 



- ) = \h J *f"^o) + \s 2 J{l - Z c 2 )[25 mn - 5 m+ln - 5 mn+1 ] . (5.32) 



Obtaining the actual numerical coefficients in eqs. (??) and (??) is somewhat tedious; 
however, the fact that the m, n-dependence enters only via 5-functions of the form S mn+ k, 
which is all that is needed for the following argument, can be anticipated directly from the 
fact that the original C of eqs. (??) to (??) is "translationally invariant" (i.e. invariant under 
I — > I + 1). Now, when writing 5 2 C of eq. (??) in terms of the Fourier transform variables 
(5&i,5Zi) of eqs. (??) and (??), one may utilize the relation 

5 5 

£ 5 mn+k 5z m 5z n = e^e'^^, , (5.33) 

mn=0 1=0 

which is essentially a consequence of the convolution theorem for discrete Fourier transfor- 
mations. Using eqs. (??) to (??) in eq. (??), it follows directly that the variables (5$ ,dZ ) 
decouple completely from the others in 5 2 C [?]. Hence the path integral (??) factorizes into 
a product, say Iih, of two independent path integrals. l\ depends only on (8& ,8Z ), the 
variations of the tunneling coordinates around their tunneling path values: 

J 1 = Jv(6n o )e-f dTS2 £°/ h . (5.34) 
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In writing this equation, another shortcut has been employed: 5 2 C is the second variation 
of the simple Lagrangian C Q of eq. (??), namely 



S 2 £ = -ms5$ 5Z + 12s 2 J5Z 2 + \ J b f '($„) 5$ 2 . 



(5.35) 



This simple form may be used since 
Ansatz (??), is identically equal to 



9 2 £ 

d*Z 



, evaluated according to the tunneling path 
with similar comments holding for the other 



\d 2 Z 

second derivatives of C Q with respect to the $ and Z . 

The other path integral, J 2 , depends only on (8&i, 5Zi)i^ , the variations of the transverse 
coordinates around their zero values along the tunneling path. These fluctuations are exactly 
harmonic oscillators in the limit -jij — > 0, since then all time dependence drops out of the 

coefficients in eqs. (??) to (??) (recall that Z c = iZ Q + z g , with i5Z D oc \Jb/a oc \J Jb/s 2 J 
). More generally, the fluctuations of the transverse coordinates will shift the ground state 
energy by an amount proportional to their characteristic frequencies. However, they will 
not make a contribution to the tunnel splitting energy A. 

Thus, for the calculation of A, only the path integral Ii is relevant. Since this depends 
only on C , the demonstration that only C Q is important for the calculation of A is complete. 

It follows that A may be directly obtained from the effective Lagrangian L Q of section ??, 
eq. (??), inserted into the results of section ??. Substituting the values of eq. (??) and (??) 
into eq. (??) gives 







7T 2ip a 



;j 6 /"(0))ij2exp [-^/3J b /Jf^f(<p) 
if s = half-odd-integer . 



if s = integer , 



(5.36) 



This is the main result of this paper - an explicit expression for the tunnel splitting 
energy in a non-trivial setting. 

Expression (??) can be evaluated numerically (see the last paragraph of appendix ?? for 
the definition of J b and f(tp), and eqs. (??), (??) and (??) for <^ as ). The results for A, in 
units of J, for some values of s are tabulated in table ??. It should be remarked, though, 
that strictly speaking and for several reasons, our methods are only applicable in the limit 
of s ^> 1. Collecting the s-dependencies of y2 as , /"(0) and J b gives 



A oc sJexp[5.6s 1/6 - 5.6s 1/2 ] for s > 1 



(5.37) 



s 


1 


2 


3 


4 


5 


6 


10 


20 


A/J 


2.3 


0.83 


0.32 


0.13 


0.06 


0.02 


0.001 


5 xlO" 6 



VI. DESTRUCTIVE INTERFERENCE 

It is quite striking that the tunnel splitting energy and the tunneling amplitude (see 
eqs. (??) and (??)) are exactly zero if s is a half-odd-integer. This is a consequence of 
destructive interference of the (+)-instanton and (— )-instanton amplitudes in eq. (??). Ac- 
tually, this result can be arrived at without the need for a detailed calculation (see also [?]). 
All that is needed is the fact that Re[S^] = Re[S~] and Af + = Af~, which is guaranteed by 
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symmetry (under some technical assumptions, explained in section ??), and knowledge of 
the phase of the classical action, eq. (??). In fact, a similar result holds for a much larger 
class of tunneling problems on the Kagome lattice. 

Consider the simultaneous tunneling of larger sets of spins on the Kagome lattice. Take, 
for example, any closed "loop" of spins such that all the spins on the loop alternate between 
the type A and type B around the loop. It is proven in appendix ?? that all such loops 
contain 4n + 2 spins (with n some integer). Now consider a "generalized weathervane mode", 
i.e. tunneling between configurations \i) and |/) that differ from each other only through 
4>i —* 4>i + 7T f° r eacn spin on the loop (i.e. in that all A- and B-spins on the loop are 
interchanged). This generalized weathervane mode does not cost any classical exchange 
energy, but selection effects provide a coplanarity barrier that has to be tunneled through. 
Under the (admittedly somewhat tenuous) assumption that during the tunneling motion 
of the 4n + 2 spins all other spins will remain fixed, the following assertions can be made: 
The absolute values of the tunneling amplitudes of the (+)-instanton and (— )-instanton 
events are equal, and the relative phase between them is i2its(An + 2)(— |). Hence, for half- 
odd-integer s, destructive interference between (+)-instantons and (— )-instantons occurs 
again. 

Only when two loops, with a total number of 4n spins, tunnel simultaneously and syn- 
chronously will there be constructive interference between all the 4n spins. The smallest 
group of spins for which this is conceivable is the 6-spin inner loop of a hexagon, nested 
inside a larger loop of 18 spins. These two loops would tend to tunnel synchronously, be- 
cause of the coplanarity forces between the triangles connecting them. Thus, the smallest 
tunneling event that has a non-zero amplitude involves a double loop of altogether 24 spins. 

The fact that destructive interference happens for half-odd-integer s is noteworthy, since 
s is half-odd-integer for the two experimental realizations of the Kagome lattice that have 
been proposed. These are the magnetoplumbite SrCrs- x Ga4+ x OiQ (s = §) [?,?,?] and the 
second layer of He 3 atoms on a graphite substrate (s = \) [?,?]. 

VII. CONCLUSION 
A. Discussion: Comparison of energy scales 

In this paper, the tunneling barrier was taken as a fixed potential, given by H cop in 
eq. (??), although the original Hamiltonian (??) had no "coplanarity forces". In what 
follows, we argue that this approach can be justified in the limit of large s. 

Our estimate of H, cop in appendix ?? is based on the following strategy: The tunneling 
hexagon is frozen at a given point along its tunneling path, at which its plane is tilted by 
tp relative to the reference plane. Then the frequencies u>i(tp) of the zero-point modes of the 
spins on neighboring hexagons are calculated as a function of tp. The coplanarity barrier is 
then proportional to J2i iM^Kv 9 ) — ^z(O)). 

This procedure is reasonable if the tunneling process is "adiabatically" slow, i.e. if 
the duration of an instanton, say r in , is indeed long compared to the typical oscillation 
periods cu l (tp) of the zero-point modes on the surrounding hexagons. In this case, we are 
effectively "integrating out" in a crude way the fast modes of a complicated many-spin 
problem. However, if for some zero-point modes one has uj 1 (tp) ~ r in , these modes are 
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slower than the instanton and do not have time to affect the barrier. In that case, the 
effective barrier should be smaller than the one we used, since the contributions of slow 
modes to TL cop should not be counted. 

Now, in appendix ??, the behavior of the zero-point modes of the hexagons surrounding 
the tunneling hexagon is found to be as follows: There are both ordinary and soft modes 
(corresponding to pf^ 3 and pf =3 , respectively, in appendix ??). For the ordinary modes, we 
find hui{ip) ~ sJ, as is usual for antiferromagnets. The soft modes are weathervane-like 
fluctuations; as tp ~ 0, we find /k^ so ft(y?) ~ s$J, while for large </?, one has hu) so ft(<p) ~ sJ. 
If one were to consider fluctuations on the entire lattice (instead of only on the neighboring 
hexagons of the tunneling hexagon), these soft modes would translate to an entire brach of 

2 

soft modes, for which one would also expect fta> so ft ~ s^J. This, indeed, is the result found 
by Chubukov [?], using a method based in Fourier space. 

The duration of an instanton event, on the other hand, scales with s according to r; n ~ 
1/d ~ Ti/(Jsi), according to eqs. (??) and (??). Thus, in the limit of s >> 1, one has 
T in uj 1 for all zero-point modes, both ordinary and soft. In this limit, the instanton is 
thus indeed slow compared to the zero-point fluctuations, justifying the use of a coplanarity 
barrier H cop for sufficiently large s. 

How does the calculated tunnel splitting energy A compare to the in-plane selection 
energy (of order J c ) mentioned in section ??, that tends to lift the degeneracy between \i) 
and |/), and hence to suppress tunneling? According to Chubukov [?] (see also [?]), J c is 
a small fraction of the zero-point soft-mode energy, i.e. J c ~ r/Js 2 ^ 3 , with rj < 0.05. Since 
we found = 5/2Js (see eq. (??)), neglecting J c relative to Jf, is justifiable, particularly 
when s ^> 1. This is the justification for not adding an in-plane selection term to 7ihex i n 
eq. (??) in the same way that we added the coplanarity selection term H cop . However, it 
should be noted that for small tunneling angles if, for which Ti cop ~ Js 2//3 , the inclusion of 
in-plane selection effects, also of order Js 2 / 3 , would probably affect some of our results. For 
example, our calculation of y? as in section ?? would be affected if Jsip* < J c < 

Finally, from the table on p.??, it follows that for physically realizable (i.e. small) spin 
values, we have A ~ J c . Thus, to the extent that our calculations can be trusted for 
small s, tunneling effects are at least as important as in-plane selection in determining 
the true nature of the ground state. We conclude, therefore, that for small (integer) s, 
tunneling should be regarded as a significant disordering mechanism competing with order- 
from-disorder selection effects. 

In a quantum Potts model of the entire Kagome system, the Hamiltonian would have 
matrix elements proportional to A between any two pairs of states which differed by the 
exchange A <-> B along a single 6-spin hexagon loop (or more generally, there would be 
matrix elements for any larger closed loop of alternating AB-spins). Then a superposition of 
Potts states, which takes advantage of the resonance between different hexagons (or larger 
loops), would have a lower energy (at small s) than the a/3 x \f?> state which is favored by 
the J c term. 

Two possible kinds of ground states suggest themselves: The first is a "totally disordered" 
Potts state; it can be idealized by a trial wavefunction which is an equal superposition of all 
the Potts groundstates, which is known to be disordered (it has only algebraically decaying 
correlations [?]). The second kind of ground state is a "partially disordered" a/3 x a/3 state; 
an idealized trial wavefunction has, say, all C-spins fixed, as in the usual \/3 x \/3 state, 
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but each ABABAB-hexagon is replaced by an equal superposition of the ABABAB- and 
BABABA-configurations, (with no correlation between one hexagon and another.) Such a 
state would still have long-range order with respect to the C-spins, but would be disordered 
with respect to the A- and B-spins. We call this a "ferrimagnetic" ground state, since the 
expectation value of n z is +1 on the C-sublattice, and — \ on the other two sublattices. 

Both these candidate states have considerably less order than the long-range-ordered 
\/3 x \/3 state. On the other hand, they are more ordered than the "spin-nematic" state 
which has been previously proposed for this system [?,?]). In the "spin nematic", all spins 
lie in the same plane, but all directions in a given plane are equally likely, i.e. (e t9 ^) = 0. 
The states described by quantum Potts models, have either an A-, B- or C-spin at each site, 
so that (e m ^) ^ [?]. 

For half-odd-integer s, of course, hexagon tunneling is absent (the smallest tunneling 
event with non-zero amplitude involves the synchronous turning of 24 spins), as discussed 
in section ??. Therefore, one would expect the long range order of the usual a/3 x a/3 state 
to persist, at least to much smaller values of s than for integer spins. 

B. Possible extension and other applications 

A more fundamental approach for treating coplanarity selection effects, that should be 
trustworthy also for small s, would be as follows: write down a functional integral for 
the entire iV-spin system, starting from the Heisenberg Hamiltonian, eq. (??), without an 
explicit coplanarity term such as eq. (??). Adopt a tunneling path in which only the six 
spins of the tunneling hexagon move. The tunneling coordinates would be essentially the 
same as the collective coordinates ($ D , Z ) introduced in section ??. However, one would 
then have of order N — 1 transverse modes, instead of only 5. Integrating these out would 
result in an effective potential for (f ,Z ), which could, in principle, be non-local in time 
(since it encodes the effect of changing ($ G , Z Q ) at time r on the transverse modes, which, 
if these modes are slow, will in turn influence the potential felt by (<3> , Z a ) at later times). 

We briefly mention a possible application of tunneling calculations in the context of 
finite-size systems. Consider some frustrated antiferromagnetic system whose ground states 
have spin order with non-trivial discrete symmetry. For small system sizes, the eigenvalue 
spectrum can be found by exact diagonalization of the Hamiltonian. Tunneling calculations 
could be useful in obtaining better interpretations of such an eigenvalue spectrum. For exam- 
ple, the Ji — J2 square antiferromagnet (for J 2 > 0.5 Ji) has classical ground states in which 
the even sublattice has alternating spins of one staggered magnetization and the odd sub- 
lattice has an independent, arbitrary staggered magnetization. When quantum fluctuations 
are taken into account, spin- wave theory gives an effective "collinearity" force (analogous to 
the coplanarity term in eq. (??) in the present paper). This lines up the staggered magne- 
tizations along the same axis, in either sense [?,?], so that classically there are two ground 
states related by a discrete symmetry and separated by a barrier. Consequently, the two 
lowest eigenstates of a large system should both be spin singlets, namely the symmetric and 
antisymmetric combinations of the two classical ground states (which are also spin singlets) 
that have been mixed by quantum tunneling. (This behavior contrasts to that of ordinary 
antiferromagnets; there the smallest gap is associated with E ~ where y is the sus- 

ceptibility, so that the first excited state is a triplet.) The classical path expected (in a small 
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